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1 Introduction 

Duality plays an important role in the analysis of quantum field theory (QFT) and statistical 
physics. This property allows to study the behavior of observable in the strong coupling region 
of one model in terms of the ones in the weak coupling region of the other (dual) model. In 
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different region of the coupling constant, it is then possible to use perturbative and semi classical 
methods. For instance, in four dimensions the electro-magnetic duality conjectured in and 
developed in p[ is the main ingredient in studying the spectrum and phase structure in = 2 
supersymmetric Yang-Mills theory. 

Almost all two-dimensional relativistic theories can be understood as conformal field theories 
(CFT)s - describing a fixed point of the theory - perturbed by a relevant operator. If the 
perturbation preserves integrability, then the analysis remains exact beyond this point. In such 
cases, the non-perturbative analysis essentially simplifies. Besides its Lagrangian formulation, 
the QFT also possess an unambiguous definition in terms of factorized scattering theory. These 
data permit one to use non-perturbative methods for the analysis of integrable QFTs and make 
it possible, in some cases, to justify the existence of two (dual) Lagrangian representations of the 
theory. 

Among two-dimensional dual models, one finds the well-known sine-Gordon/massive Thirring 
QFTs 0. Another interesting example of duality is the weak-strong coupling duality flow from 
the (affine) Toda theories (A)TFTs based on the (affine) Lie algebra Q to the theory based 
on the dual (affine) Lie algebra Q M. More generally, integrable deformations of ATFTs also 
provide series of dual models |Q |, ^possessing many applications. For instance, the scattering 
data (S'-matrix) for the ones based on Lie superalgebras with massive excitations only have 
been considered in detail |[. These models correspond to ATFTs (for the purely bosonic 
part) coupled with one or two Majorana fermions. This was done either for models based on 
A^'^\0, 2r — 1) and C^'^\r+1) (whose bosonic root subsystem is simply laced) but also for B^^\0, r) 

and A(^)(0,2r) (non-simply laced subsystem corresponding respectively to BCr and Ci^^ roots). 
In particular, these cases can be obtained as reductions of more general deformed ATFTs 0. 
Furthermore, they can be seen as integrable perturbations of CFTs with extended symmetry 
based on WB{0,r) (or also called fermionic WBr) algebras. 

Let us now consider some of the simplest cases above, i.e. corresponding to a rank r = 1. 
Then the models C^'^\2), B^^\0, 1) and A^'^\0, 2) can be understood as integrable perturbations of 
WB{0, 1) CFT, i.e. the = 1 superLiouville (SL) field theory. Integrability can be shown either 
by using the so-called singular vector analysis |^ or by constructing explicitly local conserved 
quantities (see for instance for the two latter cases). In the first model C*-^^(2) the perturbation 
preserves supersymmetry and it corresponds to the well-known supersinh-Gordon model. In the 
two other cases, which will be the subjects of the present analysis, supersymmetry is broken. 

If the theory contains only massive particles (like those presented above), the S'-matrix data 
exhibit information about the long distance (IR) property of the theory. On the other hand, 
the knowledge of the CFT data and the relevant operator associated with the perturbation 
also define completely the theory. Once identified, the CFT contains all information about the 
UV behavior of the theory. It is consequently important to connect these two kinds of data 
in order to understand better the structure of integrable QFT. When the basic CFT admits 
a representation of the primary operators in terms of vertex operators the CFT data contain 
the so-called "reflection amplitudes" , which relate different vertex operators possessing the same 
quantum numbers. Considering the system on a circle of size R, it has been shown fl^ that 
reflection amplitudes play a crucial role in the determination of the effective central charge 
Ceff{R) in the UV domain {R — * 0) of the QFT. On the other hand, the effective central charge 
Ceff{R) can be calculated independently from the S'-matrix data using the TBA method. At 
small R, asymptotic can be compared with the one following from CFT data, if one knows 
explicitly the exact relation between the mass of the particles (IR data) and the (UV) parameter 
in the Lagrangian. The agreement of both approaches can be considered as a non trivial test 
for the S-matrix elements. On the other hand, reflection amplitudes are a key tool in the 
calculations of exact vacuum expectation values (VEVs) in integrable QFTs. It has also been 
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shown that perturbative analysis and semiclassical methods support these conjectures. Both 
types of non-perturbative analysis have been performed in many cases |jlO|, 0, |12|, [l^, |15|, [T^ . 
Beside the confirmation of the S-matrix and calculation of exact quantities, they also provide a 
useful tool to check duality. 

Perturbed CFTs have also attracted much attention since they can describe many physical 
systems in the vicinity of a critical point. Several models connected with either ATFTsQ or 
integrable deformations of ATFTs have received particular attention in the recent years, for a 
large class of strongly interacting solid state physics problems. 

In a recent paper we considered integrable coupled identical minimal models for which 

1^ for which 



the interaction preserves integrability. These kind of models were introduced in |]J 

exact results like VEVs and mass-parameter 

ATFTs 



17 



.(1) 



the scattering properties were considered. In 
relations are obtained using exact results for 

In this paper, we consider in detail the ATFTs based on super Lie algebras -B*^^''(0, 1) and 
y4*^^^(0, 2) which correspond to a critical Ising model coupled with a Liouville field theory. Using 
non-perturbative analysis based on reflection amplitudes, in the UV limit of both models the 
effective central charges are calculated in the Neveu-Schwarz (NS) sector (subsect. 2.1). Exact 
VEVs of primary fields which belong to the Liouville field theory and their first descendants are 
proposed (sect. 3). Observables of both models are shown to interchange with the flow of the 
coupling h. In particular they coincide at the self-dual point 6^ = 1. Duality relations between 
parameters in both Lagrangian are proposed as well as exact relations with the mass of the 
particles. 

In section 4, taking analytic continuation of the coupling h ^ ij3 and fixing its value, several 
kind of coupled minimal models can be obtained. Namely, a critical Ising model and minimal 
model ^A.p|pl interacting through e(x)$i2(a;) or e{x)^2i{.x) where e{x) denotes the energy operator 
of the critical Ising model and $12 (a^) and $21(2;) are primary operators of the second model. 



For instance, together with the results of |]T^ this completes the analysis for Ising-Ising coupled 
models, as energy-spin interactions are considered here. Among other models, exact results 
are obtained for the Ising-3-state Potts models coupled by energy-energy which appear in the 
phase diagram of Zg spin models [^]. Furthermore, we will relate the critical and tricritical Ising 
models coupled through this type of interaction to a special case of the homogeneous sine-Gordon 
model (HSG) SU{?,)2 



1 see 



2^ and references therein). This model which can be understood 
as an integrable perturbation of the WZNW-coset CFT S'f/(3)2/(C/(l))^ have recently attracted 
attention [^, |2^, |23|, It possesses applications in one- dimensional quantum spin systems: 
such perturbations appear directly in the CFT study of a three-leg spin ladder which consists of 
three spin-1/2 Heisenberg chains weakly coupled by on-rung interaction p5| . 



However, there are also other coupled models which, although studied numerically, have not 
yet been studied analytically in great detail. For instance, the coupled Ising-XY model has 
attracted great attention as it is expected to describe the critical behavior of a large class of 
two dimensional classical XY models having the particularity to exhibit both discrete Z2 and 
continuous U{1) degeneracy in their ground state. Among these models, one finds the fully 
frustrated XY model (FFXY) g, J1-J2 XY model pf, triangular lattice frustrated XY 
model |2^, helical XY model |3^, the Coulomb gas system of half-integer charges |^,... In 
particular, the FFXY model can be physically realized as a Josephson-junction array of large 
capacitance in a perpendicular magnetic field corresponding to a half-fiux quantum per plaquette 
1 3^ . For imaginary values of the coupling, both models considered here^ possess such Z2 x f/(l) 



^For instance, the analytic continuation of An-i ATFT provides exact results for leading thermal perturbation 
of Z„-parafermionic CFTs Also, integrable perturbations of minimal models have been considered 

■^The model associated with C*^^^(2), which corresponds to the = 1 supersymmetric sine-Gordon model, has 
been suggested as a candidate for the coupled Ising-XY model by Foda [M. Exact results for this case follow 
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symmetry, so that we apply previous results in section 5. 

The models ^^^■'(0, 1) and ^^^•'(0, 2) describe integrable perturbations of WB(0, 1) supercon- 
formal minimal models, i.e. superLiouville theory. At specific values of the coupling, one can 
then obtain VEVs and exact mass-UV parameter relation in = 1 superconformal minimal 

models perturbed by Neveu-Schwarz operators $1^"^^ ^31^^ '^15^^'* • For simplicity, only the 
first case is considered in detail in section 6. Some checks are done which support the results. 
Finally, some remarks about the Ramond (R) sector of both models follow in the last section. 



2 Liouville field theory coupled to Z^-parafermions 



Instead of studying coupled Ising-Liouville models we rather prefer to consider two different kinds 
of interactions between a Z„-parafermionic and a Liouville field theory. Although probably 
not integrable for generic values of n 7^ 1,2, it provides a useful way to check the consistency 
of the expressions for n = 1,2. Consequently, we consider two QFTs which admit Lagrangian 
representations described in terms of a parafermionic field theory interacting with a Liouville- like 
term 



A. 



in) 



-bip _i_ ~ 2b(p 



(1) 

(2) 



Here aIt^^ denotes the Z„ parafermionic CFT with central charge c 



,Q i^ciiuuco tiic ^„ paiaiciiiiiuiin. ^j. ± vviuii i^ciiuicii i.iiaig,c o — 2 — and the ficlds 1p{lp) 

are the holomorphic (antiholomorphic) parafermionic currents with spin s = 1 — ^ (s = — s). 

For n = 1, we have iplip) = I in ([^), : the Lagrangians above coincide with the (integrable) 
well-known Sinh-Gordon (ShG) and Bullough-Dodd (BD) model, respectively. For n = 2, the 
parafermionic current ijj is a. Majorana fermion and to cancel fermion divergencies we have to 
add a counterterm /ie-^^^ in both cases (|l|), (^. The resulting Toda-like part in the QFTs (|I]), 

then becomes respectively a BD or a ShG model. It is known that these QFTs are integrable 
for n = 2 for which the factorized scattering theory has been studied in detail in p. 

Like for most of two dimensional field theories, conformal perturbation theory (CPT) can be 
used to study these models. Indeed, QFT with action (|l|) (or (|) with the substitution n k) 
can be seen as two different deformations of the same following CFT 



X 



in) 
CFT 



(3) 



which, for n = 1,2, coincides with the Liouville and A^ = 1 supersymmetric Liouville models 
(adding the counterterm), respectively. The stress-energy tensor 



j^in). 



^(9^)2 -g'^V with Q' = \ + ^^^ 



(4) 



and Tp"'' associated with the parafermionic CFT, ensures the local conformal invariance of (^. 
Here, we denote d = |(9^ — idy) and the fields are normalized such that 



Gaussian 



— In |a; — y\ 



(5) 



from 1 16 so we will not discuss this model here. 
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For instance, we have the conformal dimensions 



i-i 

n 



a{2Q' + a) 



(6) 



Besides the conformal symmetry, for general values of n the parafermionic CFT possesses 
an additional symmetry generated by the parafermionic current 

in this CFT are the order parameters aj, j = 0, l,...,n — 1 with conformal dimensions 



The basic fields 
6, = 



j(n ~ j) / {2n{n + 2)). The operator algebra also contains Z„ neutral fields e^^', j 



< n/2 



with conformal dimensions Dj = j(j + l)/(n + 2). All fields in this CFT can be obtained from 



the field aj by application of the generators of the parafermionic symmetry [^]. It is thus natural 
to introduce as the basic operators the local fields aje"''^. However, here we will essentially focus 
on the Liouville part, i.e. we will take j = 0. 

The exponential fields Va = e**"^ are spinless conformal primary fields of the CFT (|^) 
conformal dimensions —a{2Q' + a)/2. It can be shown |T6[ that the fields Va and V-2Q' 
reflection images of each other, i.e. they are related by the linear transformation 



with 

a are 



= 4") (a) e-(2«'+'^)^ 



(7) 



For n = 1, the reflection amplitude Rl (a) reduces to the so-called "Liouville reflection am- 



plitude" proposed in |T^. For = 2, it corresponds to the N = 1 superLiouville reflection 
amplitude in the NS sector proposed in [Q. For general values n, the (NS for n = 2) reflection 



amplitude 



a) associated with action @ writes ||16 

2(a + Q') 



Rt'\a) = - 



7TK'-f{Q'b) 



n 



■n264x l/n 



r(l - (a + Q')h)V{l - 2{a + Q')/nh) 



r(l + (a + Q')h)V{l + 2(a + Q')/nh) 
Here we denote 7(x) = r(a;)/r(l — x) as usual. 



2.1 Scaling functions for generic values of n 

Consider now the CFT (j^) on an infinite cylinder of circumference 2tt with the cartesian coor- 
dinates Xi, X2 where X2 along the cylinder is defined as the imaginary time and xi ~ a;i + 2tx 
is the space coordinate. For simplicity, we will here focus on the sector in which only bosonic 
zero-modes appear. In particular, for n = 2 it corresponds to the NS sector. Then, the reflec- 
tion amplitude (P) provides non-perturbative information in the study of quantum mechanical 
problem for bosonic zero modes 

['^'^ dxi 

ipo = / ip{x)— (9) 

Jo ^TT 

of the fields (f{x). In the semi-classical limit 6 — 0, where one can neglect the oscillator modes 
of <f{x), the Schrodinger equation governing the zero-mode dynamics writes Q 

n^^^^p{ipo)=Eo^pM with Eo = -'^ + P^ (10) 
^For n = 2, one has for instance |l| Ti'^^ = - (af^)^ + 27r/ie-2''Vo 
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the ground state energy where the momentum P is a real vector. The full quantum effect can 
be implemented simply by introducing the exact reflection amplitudes which take into account 
also non zero-mode contributions |jTO|. 

The wave function in the asymptotic region can be obtained using the same arguments as 
the ones already applied for other models. Namely, the exponential term in the Hamiltonian is 
considered as a potential wall. An incident plane wave with momentum P is then reflected to 
the plane wave with —P. The phase change corresponding to this process is associated with the 
reflection amplitude given above. Consequently, the wave function \l/p((y9o) is simply written as 
a superposition of two plane waves: 



in), 



'iP(po 



at 



+00 . 



with Sl"'\P) = R^\—iP — Q'). Using the approach proposed in |]TOl, we can now obtain the 
scaling functions in the UV region of the QFTs (|l]) or (|^) defined on a cylinder with circumference 

Let us first consider the QFT ([l|). The additional term /xe'"^ in its action compared to the 
CFT introduces a new potential wall, which implies the quantization of the momentum P in the 
wave function. It depends on the size of the enclosed region, which is proportional to ln(l/i?). 
As we will see later, this quantized momentum P{R) defines the scaling function Ceff{R) in the 
UV region using eq. (p^Ol). 

In action (|l]) the dimensions of the parameters are Dim[K] = f + and Dim[fi] =2 + 6^. It 
is now convenient to rescale back the size of the system from R to 27t. The action (|1|) becomes 



T = l 



An) 



ipipe 



Following previous analysis (see JTO], |T2| for instance) and due to the form of the perturbing term 
in (|T|) we obtain the following quantization condition for the momentum P: 



, R ^ 



-^2iP{b+l/b)Hn 0(1) 



si,''{p)u^-,sir\p) = i 



(12) 



where Sl^\P) is nothing but the so-called "Liouville reflection amplitude" proposed in |]T0|. For 
further convenience, here we introduced the "deformed" Coxeter number ||37|| : 



Hr, 



2(n + l] 



n 



'l-B) + 2B 



with 



B 



1 + 62 



(13) 



For the lowest energy state, in terms of the reflection phases, eq. (|12D reduces to 



LP = 2ti - 5\^\P) - 5^^\P) 



(14) 



with 



|i/„(l + f=)ln(^; 



In 



7r/i7(672) 



2 ,^,7r«:7(l/n + 6V2)^2 



n 



(n?h^y/r. 



and we used the convenient notation 



5^^\P) = 



T{l+iPh)T{l + i2P/nh) 
\.V{l-iPh)V{l-i2P/nb)\ 
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In the UV region i? — we can solve eq. (14) perturbatively by expanding the reflection phases 
in powers of P. We obtain 

iP = 2n- (<5g + 5g)p3 - + 5g)P^ + ... (15) 

where we define 

and introducing the Euler constant ■je, 

£^L-Lo = L-27.^^^^i±^. (16) 

b 

The ground state energy of the system on the circle of size R is given by E{R) = —7iCeff{R)/6R 
with the effective central charge Ceff{R) = {n + l)/2 — 12P^. Here, P is the solution of the above 
quantization condition (14), and perturbatively (plSf ) in powers of After some calculations, 
we find that the UV behavior of the QFT (|l]) - in the NS sector for n = 2 - is characterized by 
the scaling function 



,27r 2 2C(3)(263 + 8(n=^ + l)/n363) 27r, 
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^ 2 C(5)(2&5 + 32(n^ + l)/n5&5) ^27r^, ^ e((^'")^) . (17) 



In particular, for n = 1 this result agrees perfectly with the sinh-Gordon one |T0| . 

The same analysis can be performed along the same line for a different kind of perturbation, 
i.e. for instance in QFT (0). In this case, the dimension of the parameters in the action are 
Dim[K\ = ^ + b'^ and Dim\p] = 2 + 46^. The quantization condition is now 

^^)-2.P(.+l/.)H„^(l)(p)|___^(n)(^^ = 1 , (18) 

where we define the "deformed" Coxeter number Hn = ~^(1 — -B) + 3B. As before, eq. ([I8|) 
can be written in terms of the reflection phases as LP = 2it — 6^fj (P) — 5^"'' (P) with 

L = -lUl + h-)\n{^) - lln [[vr/27(26^)][^ ^^'/? + ^'/'V j • (19) 

Consequently, it is straightforward to show that the scaling function for the UV behavior of the 
QFT (Ef) - in the NS sector for = 2 - is given by the following expansion 



n + 1 



,27r 2 _ 2 C(3)(9&=^ + (8 + n=^)K&=^) 27r 5 



^ 2 C(5)(33&5 + (32 + n5)/n565) ^27r^, ^ 0((^)')] (20) 

with £ similar to £ in (|T6]) but with the replacement Hn —>■ Hn- It can be checked that for n = 1, 
this result coincides exactly with the Bullough-Dodd ||T^ scaling function, as expected. 

Notice that for the specific value of the coupling constant 6^ = 1, both scaling functions take 
the same form for any values of n. 
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2.2 Thermodynamic Bethe ansatz for n = 2 and duality 

For specific values of n, the effective central charge calculated above from the CFT data (reflection 
amplitudes) can be compared with the same function determined from the numerical solution 
of the TBA equations for the QFT ([l|) and (^. For n = 1, this analysis has been done in [|in] 
for the ShG and in [0 for the BD cases. Then, let us consider the case n = 2 in the QFTs 



and (0). As was conjectured in ||^ (even for higher rank r > 1 of the affine Toda part) these 
QFTs possess a weak-strong coupling duality: the analysis of the factorized scattering theory 
shows that there exists a QFT which possesses two (dual) perturbative regimes associated with 
action (|l]) and (^, respectively. An intermediate mass spectrum (consisting of two particles for 
the rank r = 1) was proposed: 

M^ = m and Mi = 2msm{n / H2) . (21) 

This mass spectrum flows from the one of the QFT (|l]) to the one associated with the QFT @ 
while b increases. 

Using the corresponding (diagonal) FST proposed in p - we report the reader to this work 
for details - one writes the TBA equations. Namely, 



-el(e,i^)^ 



d9 (22) 



c^eff^\R) = ^ J cosh^ [m^ ln(l + e-^-^^^'^)) + Mi ln(l + e 
where the functions e^{9,R), ei{9,R) satisfy the system of 2 coupled integral equations 

M,Rcoshe = e,{e,R)+ J" [ ip.JO - e')\n{l + e"'^'^''^^)^ for 

^ — ^ / 27r 

with the kernels ipij defined as the logarithmic derivatives of the ^-matrix elements obtained in 
However, the function E^'^^^\R) defined from the TBA equations differs from the ground state 
energy E{R) of the system on the circle of size R by the bulk term E^'^^^^R) = E{R) — fi-'^^R, 
where fr^^ is a specific bulk free energy of the QFT. To compare the same functions we should 
then substract this term from the function E{R) defined by cfjj{R)^ i.e. for instance in case (|l]) 
we have 



cSr(ii) = c3',(fl) + 5^/a (23) 



and similarly for case (|^). Notice that the contribution of bulk term fr becomes quite essential 
at -R ~ In many examples of known QFTs, this quantity can be calculated explicitly using 

Bethe Ansatz (BA) method (see for instance PBl |3^). Here, we propose the following expression 
for the specific bulk free energy in the QFT ([1|) (the same quantity for the QFT (^) follows using 
the duality 5 ^ 1 - 5 and iJa ^2): 

.(2) ^ ^ sin(7r/g2) , 

^"=^ 8 sin(7r5/ff2)sin(7r(l-5)/i/2) ■ 

Finally, to compare the expansion coming from the CFT data and the one from the numerical 
analysis of TBA eqations, we need the exact relations between the UV parameters in the actions 
(P and (@) and the IR mass scale for the particles m. Here, we propose 



-vr/.7(l + &V2)f[f7(^)]^ 



mr(i + ^)r 



2l+2/H2r(l/if2) 



2/^2(1+62) 

(25) 
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for the QFT (|I]) and 



+ t)r( 
2i+2/H2r(i/#2) 



l-BM 2i?2(l+fe2) 



(26) 



for the QFT (^). One should mention that under the weak-strong couphng duahty transformation 
6 1/6, these relations exchange perfectly if the parameters in action (P and (0) satisfy the 
duality relations 



7r/i7(6V2) = vr/i7(2/6' 



62/2 



and 



TT 



,1 



-71 



262 



1 ^ 


"TTK 1 + 6^ - 




[252^^ 2 



l/fe 



(27) 



Using the first relation ( P3| ) in (14) for = 2, it is now possible to compare the expansion (|T^) 
for n = 2 to (|22|) obtained numerically using (pTf) , (23) and (|^). The good agreement supports 
the approach based on the reflection amplitudes, the exact m-n-fi (and their dual) relations given 
above, the bulk free energies (^^ of (|l]) (and its dual for (H)) as well as the S-matrix elements 
conjectured in [B. Due to the following relations 



.(2) 



{R)\t 



g(2) 



i/b 



as 



1 1/6 



the weak-strong coupling duality property between the models with action (|l]) and 
proposed in [H is indeed confirmed at the on-shell level. 



(28) 
for n = 2 



3 Vacuum expectation values of local fields 

For generic values of the parameter n 7^ 1,2, there is no reason to expect the QFTs (|l|) and (H) 
to be integrable. However, for n = 1 the ShG and BD integrable models are recovered. For n = 2 
the parafermionic current ip is a Majorana fermion. As the scaling limit of the Ising model with 
zero external magnetic field is described by a free Majorana fermion field theory, both models 
describe different interactions between a Liouville field theory and a critical Ising model 0. In 
this case, it can be shown that both models are also integrable , and that respective conserved 
quantities exchange under weak-strong coupling duality 0. 

In the previous subsection we considered the UV asymptotic of the effective central charges 
in these QFTs. The calculations were based on the reflection amplitudes, i.e. using CFT data. 
However, these functions play also a crucial role in the calculation of vacuum expectation values 



in perturbed CFT []TT|. Here, at the off- shell level, we will see that VEVs also satisfy such duality 
property. 



3.1 Expectation values of primary fields and the Ising-Liouville dual 
models 

Here, we consider the expectation values of the simplest local fields which belong to the Liouville 
field theory, i.e. 

G(")(a) = {e'''P{x))r,n with r = 1,2 . (29) 

am very grateful to P. Dorey and R. Tateo for these numerical checks. 
^Notice that the form of the interaction remains the same but the definition of the Liouville coupling constant 
changes. 
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From the discussion in the previous subsection and using once again CPT framework with (|^, 
the last term in action (|l|) or (||) can be considered as a perturbation. It is then expected that 
for any r and n the VEV (|29|) satisfy the reflection relation 



(n), 



G'W(-2Q' + a) 



% V^"-; \^^^ u,; . (30) 

Let us first study the QFT (|l|). Instead of the previous picture, it is also possible to consider this 
QFT as a perturbed Liouville field theory with parameter and coupling The stress-energy 

tensor for the conformal invariant part is then 



Tl{z) 



-{difif + Qd'^ifi with Q 



b 1 
2 + 6 



Using the CPT framework, the VEV is then expected to satisfy the reflection relation 



G 



in) 



G2i(2Q-a) 



(31) 



(32) 



Obviously, it is not possible to find a solution to these reflection equations without any strong 
analyticity assumptions. Assuming that the VEV G^^2i{0') is a meromorphic function in a, the 



"minimal" solution of the reflection relations (|30|) and ( p2D can be explicitly obtained with the 
result 



G 



(n) 



7r^7(l + &V2)]'[— 7(l/r^ + &V2)]^ 



n 



-a^+2Qa 
2H„(l + b^) 



7r/i7(l + 6V2)]^^ 



X exp 



°° dt 

T 



sinh(a6t)\l/„(t, 6, a) 



sinh(t) sinh(t62) sinh(nt6V2) sinh(i7„(l + b'^)t))\ 



(33) 



where we define for — (y + ^) < TZe{ab) < y + 1 and 

^„(t, 6, a) = 2 sinh((a6 + 2Q'b)t) sinh(Q6t) sinh(nt6V2) cosh(t6V2) 
+ sinh((a6 - 2Qb)i) sinh(nQ'6t) sinhitb'^) . 

For several other QFTs (Toda, deformed Toda, ...) for which exact VEVs were proposed explic- 
itly, it has been shown that the expectation values of the fundamental field 
{ip) for 6^0 agrees with the same quantity calculated within perturbation theory. It has also 
been shown that semiclassical analysis supports the conjectures for the VEVs. Consequently, 
here we expect the same feature to be satisfied. 

We can proceed similarly for the QFT @). This QFT can be understood as a perturbed 
Liouville field theory with parameter and coupling (/i, 26). One has to consider now the stress- 
energy tensor with a form similar to (|3lD but with the substitution 



Q 



Q = b + 



2b 



The reflection relation satisfied by the VEV immediately follows 



G 



(n) 



R 



(1), 

2b > 



-a)|, 



G 



(n) 



{2Q 



(34) 



(35) 



As before, using the reflection relations ( pOD and (|35D simultaneously and assuming similarly 
strong analytical assumptions we obtain the following conjecture 



7r2 



X exp 



-462-2 ~ 



,71 K 



/l7(l + 26^)][— 7(l/n + 6V2)r 
n 



~a^+2Qa 
2H„{l + 6^)r 



7r2 



/i7(l + 26^)]-^ 



T 



sinh(a6t)\&„(t, b, a) 



sinh(t) sinh(t62) sinh{ntb^/2) sinh(if„(l + b^)t)) 
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with 

§„(t,6,a) = 2sinh((a6 + 2Q'6)t)sinh(Q6t)sinh(rf/2)cosh(t/2) 

+ sinh((a6 - 2Qb)t) smh{nQ'bt) smh{tb^) . (36) 

Here, the integral is convergent if — (y + -) < TZe{ab) < 6^ + ^ 



For n = 1, it is straightforward to check that (^) and ( p6D agree respectively with BD and 
ShG results. For n = 2, using the duality transformations of the parameters ( p7D it is also easy 
to check that 

G?i(a)|, = G?i,(a)|v.. (37) 

Finally, using the VEVs (|33D (or similarly (|36[) ) for n = 2, the bulk free energy of (|l]) can be 
calculated as we have the relation 

I^G'^M = ■ (38) 

Using (^), it is easy to show that the result for f!^J^ which follows is in perfect agreement with 
the expression proposed in (^) (and similarly for the QFT @). 

3.2 Expectation values of the first descendant fields in the Ising- 
Liouville dual models 

By adding the correct counterterms, for n = 2 both models (|l|) and (Q) can now be understood 
as two different perturbations of the N = 1 superLiouville field theory. For instance, for r = 1 it 
is given by the action 



L(^5^ + ^|J^^P) + ^{d.ipy - Kijipe-'"' + lie~M (39) 

-ZTT OTT 

and similarly with the change k k for r = 2. Superconformal transformations in the SL theory 
are generated by the super stress-energy tensor T{z) = —S{z)/2 + 9T{z) and similarly for the 
antiholomorphic part with (^, 9) the corresponding Grassmann coordinates. Here T{z) denotes 
the usual (bosonic) holomorphic stress-energy tensor given in eq. (§) for n = 2 and S{z) is a 
spin 3/2 conserved current: 

S{z) = tpdip + 2Q'd^ (40) 

and similarly for the antiholomorphic part. With the background charge given in eq. (§) for 
n = 2, the central charge of the superLiouville model is then csl = |(1 + 8Q'^). 

The basics fields in the SL theory are the operators ajexp{aip) which belong either to the 
Neveu-Schwarz sector (NS) (for j = 0) or the Ramond sector (R) (for j = 1, also called "twisted" 
fields). These primary fields have conformal dimension 

= — ^ 2 J = 0, 1 • 

Then, using the Laurent expansion 

T(z) = and S{z) = ^rz"'~^^^ (41) 

nez reZ/2 
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and similarly for the antiholomorphic part, all the other fields can be obtained via the action of 
the Neveu-Schwarz (Ramond) algebra generators L„, Sr for nGZ, rGZ+1/2 (nGZ, rGZ). 
In general, we denote these descendent fields 

where [u] = [—ui, —up] are arbitrary strings. The descendent fields (^) and the ones obtained 
after the refiection a — (2Q' + a) possess the same quantum numbers. It is then possible to 
show that the reflection property extends to all these descendants and arguments based on CPT 
approach suggest the following reflection relation for j = 0: 

{L[n]L[m]S[r]S[s]e'"^)r = Rf\a) {L\n]L\m]S[T]S[s\e~^'^^ ^°''^'^) t and r = l,2, 

where (a) is the SL reflection amplitude in the NS sector calculated in For simplicity, 
here we only consider the VEV of the simplest descendent flelds. Using the relations ( ^OD and 
( PI) above, we have 

^-i/25_i/2e'^^ = -a^^^e'^^ . (43) 

Consequently, expectation values of these operators in the perturbed theories (|l]) or (0) forn = 2 
are expected to satisfy similar reflection equations. Let us consider the ratio 

HM = . (44) 

The model with action (|1|) and n = 2 can either be considered as a perturbed SL theory or 
a perturbed Liouville theory. As before, approach based on CPT then suggests the following 
reflection equations 

H,=i{a) = W + H,^^{-2Q'-a) for a^O, 

Qj 

H,=i{a) = Hr=^{2Q-a). 

The "minimal" solution of these reflection equations is deflned up to an overall constant. To flx 
it, it is sufficient to notice that i/^=i(— 6) can be related with the bulk free energy of the system 
as follow 

H^^,{-b){e-'nr=i = ^, (45) 

which leads to the result 

(^^e = - 2(3 + 2&T ^( 3 + 2b^ ^^^sT^^ ^^-^"^ ^^'^ 

for a 7^ 0. It is straightforward to do the same analysis for the model with action (^ for n = 2. 
It can be shown that the result for {i(jijje°''^)r=2 follows from (|46|) with the substitution 6 — 1/b. 
In Appendix A, we give a further support to (^61) . In conclusion, all previous results for n = 2 
confirm at the off-shell level the weak-strong coupling duality. 
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4 Application to integrable coupled minimal models and 
the homogeneous sine-Gordon SU{3)2 model 

For n = 2, the QFT with action (0) possesses interesting apphcations. As is well known, the 



"minimal model" Aip/p' with central charge c = 1 — 6 ^^ V can be obtained from the Liouville 



pp 



model. Consequently, adding the counterterm fie ^'"^ in action (|^), the resulting QFT can be 
identified with a minimal model interacting with a critical Ising model if we substitute 

6 — > fl —ft , k ^ ik , (47) 

and make the choice either 

= Pl= p/2p' or =p^= p'/2p with p < p' . (48) 



For both values of the coupling and using the background charge (|3^) which ensures the local 
conformal invariance of the Liouville field theory, the conformal dimension of the perturbing term 
becomes A(V'^e-^^'^) = 3/372, which is relevant for (3^ < 2/3 . 

In the following, we define {$rs} as the set of primary fields with conformal dimensions 

^p^'p — ps^'^ — (^p — p')^ 
= ^ — for 1 < r < p, 1 < s < p' and p < p' . (49) 

Using the Coulomb gas representation, they are simply related with the vertex operators of the 
Liouville field theory through the relation 

$,,(a;) =iV-/exp(zr/'^V(a:)) with = + (1 _ s)/3 , (50) 

where we have introduced the normalization factors Nrs- These numerical factors depend on the 
normalization of the primary fields. Here, they are chosen in such a way that the primary fields 
satisfy the conformal normalization condition: 

{^rs{x)<l>rs{y))cFT = 1 ■ (51) 

For further convenience, we write these coefficients Nrs = N{ri^'^) where: 

N() = \- ~ ( o.2^1-^^ ^(2/3^-2r//?)^(l/2/j^ + V/3)^(2-2/3^)^(2-l/2/?^) 1| 

^ ^/^7l JJ ^r(2-2/32 + 2?7/3)r(2-l/2/?2_ry//5)r(2/?2)r(l/2/?2)J ' ^^^^ 

Taking specific values of the coupling constant in action (^, it is now possible to obtain 
non-perturbative information for two planar systems which interact through a relevant operator 
preserving integrability. The first system - described by the free Majorana fermion part - is 
identified with a critical Ising model denoted using the standard conventions. The second 

one - obtained from the QG restriction of the Liouville field theory - is identified with a minimal 
model denoted M.p/pi. Then, the action can be written either 



A = M3/4 + Mp/p> +X J d^x e{x)^i2ix) for /3 = /3+ (53) 
or A = M3/4 + Mp/p> +X (fx e{x)^2iix) for /^^ = /^^ , (54) 
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where the parameters A and A characterize the strength of the interaction and the energy 
operator of the Ising model is defined by e = iipip. Also, both QFTs make sense if 3p < 4p' in 
( ^3|) and 3p' < 4p in (^4|), which ensure the perturbation to be relevant. 

For imaginary values of the coupling h = ip, it is expected that the QFT (|^) possesses 
complex soliton solutions which interpolate between the degenerate vacua. It can be shown that 
this model possesses the QG symmetry associated with Uq{B{0, 1)^^-*) where q is the deformation 
parameter. Naturally, there are good reasons to believe that the S'-matrix of this model can 
be expressed in terms of the i?-matrix associated with this deformed affine Lie superalgebra. 
Then, in the following we assume that a breather-particle identification holds by comparing the 
resulting S'-matrix elements of the lowest-breathers (breathers with lowest mass) with the S- 
matrix elements for the quantum particles in the real coupling case. Denoting M as the mass of 
the lowest kink, we suggest the identification 



2M sin I 



2-r 



with 



(55) 



In the following, we will study successively the models with action ( p3D and (|54D . However, 
the vacuum structure in both cases is not clearly understood so the prefactor associated with 
the QG restriction is ommited for simplicity. Consequently, we denote all expectation values 
(...) = (0|...|0) where |0) is one of the many ground states. 

Let us first consider the coupled minimal models with action (|53|) . With the previous iden- 
tification and using the exact m-K-jl relations (|26|) proposed in the previous section, it is now 
straightforward to obtain the following exact M-A relations 



1-^ 



A^ 



1 lKt^)lKT^f 



1 



•2(1+0' 



7rMr(i±f) 



- ■ 



r 



■ 1 > 



2(2-S) 



for 



V 



2p — p 



(56) 



Consequently, according to ([56| ) and 
massive spectrum for: 



< 2/3, the coupled minimal models ( p3D develop a 



(^^ 



1/3<^<1, 
< ^ < 1/3, 



3m{X) 
me{X) 







I.e. 



I.e. 



1 p 

- < — < 1 

2 p' 

p 1 
< ^ < - 



(57) 



p' 



In particular, the condition (i) is always satisfied for the coupled unitary minimal models defined 
by (|53|) with p' = p+1 . Notice that for p = 2, p' = 3 the model (|53D corresponds to an off-critical 
Ising model as we have $12 = I in this case. If we make the identification fn = 2Msin(7r/3) = 
Msg where Msg is the SG soliton mass, then the M-A relation given above[] coincides with the 
one associated with the off-critical Ising model, as expected. 

Although the model (H) with imaginary coupling is very different from the real coupling 
case in its physical content (the model contains solitons and excited solitons), there are good 
reasons to believe that the expectation values obtained in the real coupling case provide also the 
expectation values for imaginary coupling. Similarly to the analysis done for the ShG and BD 



models |Tll], we obtain the VEV of primary operators which belong to the second minimal model 
using eqs. ( |36D for 6 = (|50D and ([5^ ) 



($r.(x)) 



vr^AV^ 



4-2g 



•2(1+0 



'Ml 



^ 1^ 1+g ' ' ^ 1+C 
^The parameter A is real and we obtain Msg — 27r|A| 



(1+0 . 
2-€ 



exp Q,i2((l + ^)r-2es) 



(5J 
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The function Q^uiO) for |^| < 2^ and i < ^ < 1 is given by the integral 



Qel2(^) 



t Vsinh((l + ^)t) sinh(tO sinh((2 - ^)t) 



~2t 



with 



sinh(t) 



cosh(^t) - cosh((l - ^)t) 



- sinh((2 - ^)t) sinh((l - ^)t) - cosh(t) sinh( 



cosh(^t) cosh((2 -0^) 
1 



and defined by analytic continuation outside this domain. Notice that eq. 

{<^rs{x)) = {%-rp'-s{x)) • 



satisfies 



(59) 



Finally, the expectation value ( |5^ ) can be used to derive the bulk free energy f^i2 = — lim y^^^ 

where V is the volume of the 2D space and Z is the singular part of the partition function associ- 
ated with action (|53D . The result for the bulk free energy follows from the analytic continuation 
of (HD and eq. (|55|) , i.e. 



M2 sin(^) sin(5^ 



el2 



sm 



■ 7r(l+g) - 



If we look at the action (|53[) , we can now use the relation dxfei2 
expectation value 



(60) 

(e$i2) to deduce the 



(e(x)$i2(x)) = -J 



AW(2(i^) ' 

7(S^)7(fT|: 



2-C 



1 + e 



r^(2^: 



sm 



2-^P(^)n(i5f)sin(^)sin(^) ' 



(61) 



2-5^ 



2-« 



Instead of taking $12 inside the expectation value written above, one can now consider the one 
associated with more general primary operator of the second minimal model. The corresponding 
VEV follows from the dual result of eq. (^6]), i.e. we obtain 



1+?^ J 



2(2-5) 



1 + e 



7( 



7, 



l-C+(l+g)r-2gs N 
4-2g 



2-r^2-r 



7 



- 3-g+(l+g)r-2gs N 
4-2g 



(<l>,,(x)) . (62) 



Let us now turn to the coupled minimal models with action (0). In this case, the condition 
4p > 3p' guarantees that the perturbing operator is relevant. Then, the vacuum structure is 
expected to be similar to that of (|53D . Using the substitutions 



p ^ p 



1 + e 



3^-1 

in the previous expressions, the breather-particle relation is now given by 

m = 2Msin(^^^i^) 
^ 5e - 3 ^ 



(63) 



(64) 
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and the exact relation between the mass of the hghtest kink M and A follows 



7r2 ^2rM^- 



8? 



55-3' 



r 



3g-l- 



ri 



5^-3 
2? 



(65) 



5^-3^ ^5^-3' 

the massive phase corresponds to the 



Then, for the coupled minimal models defined by 
domain: 

3 3 
(m) - < ^ < 1, 3m(X) = i.e. -<-<!. 
5 4 p' 

From eq. (|58|) and (|63D, we obtain the following expression for the VEV 



(66) 



($rs(a;)) 



5e-3 



^^AV(^)(20^ 



4€ A 



2-7(|)7(¥) 

The function Qe2i(^) for 16*1 < 2^ is given by the integral 

Qe2l{9) -- 



exp Q,2i((l + 0^-2^5) 



(67) 



t Vsinh((l + Ot) sinh(4tO sinh((5^ - 3)t) 



with 



^e2l(^,t) 



sinh((3^ - l)t) 



cosh(2^t) - cosh((2 - 2^^) + cosh(20t) cosh((5^ - 3)t) 



sinh((2 - 2^)t) sinh((5^ - 3)t) - cosh((3^ - l)t) sinh((^ - l)t) cosh((20t) 



and is defined by analytic continuation outside this domain. As was considered above, the exact 
bulk free energy can be also calculated from the results of the previous section. 



e21 



M2sin(^)sin(f±|); 



sm 



5^-3' 



(6J 



Finally, this latter expression provides us the VEV 



(e(x)$2i(x)) 



1 
A 



AW(^) 
, 7(|)7(^) 



4S 

5e-3 



4e 



^ V5?-3> 



sm 



■ 7r(3g-l) - 
• 5^-3 ' 



5e-3p(|zJ)P(^)sin(Z^)sin 



5^-3^ 



5^-3 



■ 7r(4g) ^ 
55-3 ' 



(69) 



whereas for a more general primary operator of the second minimal model, we now obtain 



(e(x)$„(x)) 



7(|)7(^; 



2? 



2g 
55-3 



^ 7(^) "L_. • (70) 



7 



■ 4g-2-(l+g)^+2g;i ^ 
5$-3 



Some checks of these expressions are desirable. For instance, the case p = 5 and = 6 in 
action (|^) corresponds to a critical Ising model coupled to a critical Za-Potts model through 
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the interaction eE where E = $21 is the leading energy operator of the Zs-Potts model. It 
is rather interesting to recall that the decoupled critical Ising-3-state Potts models appear in 
the phase diagram of the Zg spin model in the vicinity of a renormalization group fixed point. 
Using continuum field theory approach, which is valid in the scaling region around this point, 
Zamolodchikov's counting argument EOl can be used to show that the eE perturbation preserves 



conserved charges of spin ±3, ±5. As was suggested in |]20[ and confirmed here by explicitly 
constructing the QFT, this perturbation of the Ising-3-state Potts model is therefore integrable. 
Although the exact S'-matrix for these coupled models is not known explicitly, a TBA system 
based on Ei Dynkin diagram has already been considered as a good candidate The exact 



relation ( |65D as well as the breather-particle relation can been checked using this TBA system. 
The good agreement ^ supports the exact results above. 

Let us now explain how the above results can be relevant in the study of the S'[/(3)2-HSG 
model. The simplest HSG model is the complex sine-Gordon model associated with an integrable 
perturbation of the WZNW-coset model SU{2)k/U{l), whereas more complicated HSG theories 
can be viewed as interacting copies of complex sine-Gordon theories. At classical level, the 
corresponding equations of motion correspond to non-abelian Toda theories which are known 
to be integrable. At quantum level it has also been shown that integrability is preserved El 



and assuming factorization of scattering process 5- matrices have been proposed Among 
the generalizations of the complex sine-Gordon model, the S'f/(3)2-HSG can describe a WZNW- 
coset model SU{3)2/ {U (1))^ with central charge c = 6/5 perturbed by an operator with conformal 
dimension Apert = 3/5, as was shown in ||2T|. On the other hand, the case p = 4 and p' = 5 
in action (|53|) leads to a critical Ising model coupled to a tricritical Ising model (with a total 
central charge c = 1/2 + 7/10 = 6/5 in the UV limit). It is then interesting to notice that 
conformal dimensions A(e(a;)$i2(x)) = A(<l>i3(a;)) = 3/5. Obviously, in the deep UV the operator 
content of the coupled models (|53| ) is bigger than the operator content of the S'f/(3)2-HSG, 
which only consists of operators of conformal dimension 0, 1/2, 1/10 and 3/5 [^]. However, the 
operator content of (^) possesses a closed subset of operators {1, e, $12, £$125 '^'13} with conformal 
dimensions given above. Consequently, using the notations of |^ for the 5'f/(3)2-HSG, we 
propose the identification 0^'\ = $12, the trace of the energy momentum tensor 6 G {e$i2U$i3} 
whereas the remaining operator is the Ising energy density e. Then, the action associated with 
the S'[/(3)2-HSG model can be seen as a subsector of the following (integrable [^) action: 

ShSG[SU(3)2] ~ -^3/4 + M4/5 + \ j (fx^pert{,x) (71) 

with $pert(a;) = e(a;)$i2(x) + p(^i^{x) and p is a c— number. At the special value p = 0, one 
recovers (|53D . From the above analysis (|57D , it follows that the S'f/(3)2-HSG is a massive theory 
in agreement with the results of [^. Also, taking the value ^ = 2/3 in (|55|) gives m = 2M, i.e. 
the formation of stable particles via fusing of the soliton mass M is not possible in agreement 



with [^. At p = 0, its bulk free energy follows from ( |60D whereas the exact relation between 
the soliton mass M and the parameter in the Lagrangian A is given in (^) for ^ = 2/3. 

Accepting the conjectures ([58| ) and (|67D , one can then make interesting predictions for nu- 
merical values of VEVs. We report the reader to the Appendix B where various examples are 
considered. 



^P. Dorey and R. Tateo, private communication 
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5 Relation to the coupled Ising-XY model 



Over the years, the critical behavior of the two dimensional Ising-XY model, consisting of Ising 
and XY models coupled through their energy densities, has been studied numerically in some 
detail. A Hamiltonian for this model has been proposed which writes 



U/kT =-^[{A + Ba,a,j) cos(fti - e,j) + Ca^aj] 



(72) 



where A, B and C are effective couplings. The model with A ^ B is relevant for the anisotropic 
frustrated XY model [27] and anti-ferromagnetic RSOS model whereas the subspace A = B 
is relevant for the isotropic FFXY model or its one dimensional quantum version |^3[. In this 
latter case, a phase diagram has been proposed (see figure 1). It consists of three branches 
joining at a point P in the ferromagnetic region A > 0, A + C > 0. One of these branches 
corresponds to a single transition with simultaneous loss of Ising and XY order and the other 
two to separate Kosterlitz-Thouless (KT) and Ising transitions. Monte Carlo transfer-matrix 
methods |Q yields that the central charge seems to vary continuously from c ~ 1.5 near P to 
c ~ 2 at T, which contradicts the hypothesis that the line PT can be simply described in terms 
of a superposition of critical Ising and XY models with c = 1.5 as was suggested by Foda in ||3^ . 
The only possibility would be the existence of a parameter changing along the line PT that does 
not affect the symmetry. As was argued in there are three possible explanations: (a) The 
system is not conformally invariant; (6) The result is an artifact of limited strip widths; (c) It 
is a new effect. Consequently, due to the limited strip widths of the numerical analysis, there 
are some reasons to believe that the phase diagram analysis is not yet complete, in particular, 
first order transitions may appear along the line denoted PT. More recent numerical analysis of 
the Ising-XY model based on the coupled Ising-RSOS model also supports this hypothesis. 
Then, if we opt for scenario (a), it is natural to consider all possible integrable perturbations of 
a superposition of critical Ising and XY model as a starting point to study the vicinity of the 
point P. It can be shown that there exists only three kinds of integrable perturbations which 
can provide interesting candidates]^: the supersymmetric sine-Gordon model |^| which has been 
suggested by Foda or the models with action (0) and (0) for n = 2 and imaginary coupling. 



XY 










ORDERED 


DISORDERED 









Figure 1: The phase diagram for the coupled Ising-XY models 



^Notice that the continuum hmit of the generahzed Coulomb-gas representation of the FFXY model containing 
fractional charges is nothing but the (conformally invariant) action of a free Majorana fermion and a free boson. 
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Consequently, it would be rather interesting to obtain non-perturbative results for these cases. 
As some exact off-shell results (VEVS, bulk free energy,...) for the supersymmetric SG model 
16|| , here we naturally focus on action (H) for = 2|3. Using the analytic 
i(3 and /i —fl, we will consider in the vicinity of the point P of the phase 



can be found m 
continuation b — 



diagram the following action: 



As 



d X 



1 



— {ijjdijj + ijjdip) 
2tt 



Stt 



{d^ip f - k^^e-'f^"^ - 2/icos(2/5(^) 



(73) 



Before going further, let us recall known results for the XY model - which also corresponds 
to the nonlinear 0(2) cr-model. Kosterlitz and Thouless showed that spin configurations are 
mixture of topologically trivial configurations (called spin waves) and a gas of vortices with 
integer topological charge. Both are decoupled and the vortices interact through a logarithmic 
potential which is therefore identical to a two dimensional Coulomb gas. At a specific finite critical 
coupling, it has been shown rigorously [0 that this Coulomb gas possesses a phase transition. 
In the vicinity of that point, only vortices of topological charges ±1 are important and higher 
charge vortices can be neglected. In the Coulomb gas formalism, various configurations are then 



associated with the operators Oe.m with dimension de,m 



e /R + mR /A and spin s 



e,m 



em. 



In particular, the "electric" e and "magnetic" m charges exchange each other under the weak- 



strong coupling duahty R ^ 2/R. Also, the 
vertex operators in the following way 



"electric" operators can be written in terms of 



Af~^{eP)exp{ieip/R){x) with eeZ 



(74) 



where we choose the normalization factor (the mass scale M has been introduce in (^ 



7rMT{ 



2-3/32 



2 2-3^2 Yl 



1- /32 

2- 3/32 



)r( 



2/3^-1 



-1 -21 

/3 



(7r/i2^^'-27(l-2/52))i/2 



2-3/32 



(75) 



Returning to action (|7^) it is now quite natural to take P = l/\/2 in order to consider the 
Ising-XY model. For this value and R = v/2, the operators C±2,o become marginal as well as the 
last part of the action (cosine term) (^). For e > 2, operators are irrelevant. Then, we obtain 
the following action 



M 



3/4 



+ ^ / d'xid.ipY + Ap^y^ I d'x e(a;)0_i,o(a;) 



(76) 



For general values of f3, the exact relation between the mass scale M and the parameter A^ 
follows from section 2 and eq. (|55D : 



A 



7f-f\ 



1-/32 



-1 1-/3^ 



2(1-/32) 



1-/32 



/32 



(77) 



^2-3/32 /-^ V 2-3/32 . 

-'^^We have seen in the previous section that the model (||) for n — 2 can describe a critical RSOS model coupled 
to a critical Ising model through their energy densities. Furthermore, since for n = 2 both models with actions 
(|l|) and (||) possess the same limit at = 1/2, it is sufficient to focus on one only. 
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Also, using the conventions defined above, an exact expression for tlie VEV of the "electric" 
operators can be obtained: 



{Oefl{x)) 



TiMT{ 



2-3/32 



2(1-/3^) 

2 2-3/3^ ri 



X exp 



2-3/32 . 
~t 



e2/32 



'J- 12-3/32^ 



sinh(t) sinh(t/52) sinh((2 - 3/32)t) 



with 



^/3(t) 



2sinh((l - (e + l)/3^)t) sinh((l/2 - cosh(t/2) 
- sinh((l + (e - 2)/?2)t) sinh((l - p'^)t) . (78) 

For the choice (3 = 1/V2, the model ( |75D becomes massive if Xm(A^^^y^) = 0, corresponding 

to a first order phase transition. One can then check that (C±i,o) = (C'o,±i) and (C±2,o) = 
as expected. As the cosine term in the action is marginal, one can show that exactly the same 
results can be obtained if we had started from action (|l|) instead. This is not surprising as both 
models possess the same limit at = 1/2. 



6 Neveu-Schwarz sector of perturbed = 1 supersym- 
metric unitary minimal models 

Instead of considering a restriction of action (^), it is also interesting to study action (|I|) for 
n = 2 at specific imaginary values of the coupling b. As it is known, the minimal series of 
superconformal unitary models can be described from the superLiouville field theory using the 
analytic continuation 6 — i/3 in and taking the specific value 



K 



with 



K>2 



K + 2 

Their corresponding central charge is csusy = |(1 ~ k{k+2) ^ finite number of primary 

operators belonging to the NS sector are labelled by the conformal dimensions (for r — s G 2Z) 

M^s)= {r{K + 2)-sKf-A i<r<K, l<s<K + 2. 

8K{K + 2) - ' - 

Using the vertex operator representation, they can be written in terms of the exponential fields 
of the superLiouville theory as follows: 

a-r) (1 



<^l^^\x) =Af-,^exp{ir]'-'^{x)) with 



2P 



-(3 and r - s e 2Z (79) 



where we have introduced the normalization factors Mrs- Choosing a condition similar to (|5lD, 
they can be expressed in terms of the reflection amplitude in the NS sector of the superLiouville 
field theory |3§ as Mrs = M^^^Kvrs) where 



X 



2(3^'' 2 
T(l/2 + /3V2 + r//3)r(l/2 



1/2/32 - r]/p)T{3/2 - /372)r(3/2 - l/2p^ 



r(3/2 - /32/2 - 77/3)1(3/2 - 1/2/32 + r]/(3)T{l/2 + /32/2)r(l/2 + 1/2/3^) J 
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Adding the perturbing term with conformal dimension Apert 
9|), the analytic continuation b ^ ijS then gives perturbed = 



= Asiie''^) = -fe^ - 1/2 in 
1 supersymmetric minimal 



model 1^. Notice that for any values of K, the perturbation is relevant and is identified to $^^3 



(NS) 



with conformal dimension A 



(NS) 
13 



K-2 
2iK+2) 



The resulting action writes 



A 



d X $ 



(NS) 
13 ' 



X) 



(80) 



The exact relation between the parameter A and the mass scale m introduced in eq. (|21|) can be 
obtained using A = -fiAf^^^\f3) with the resuh 



A 







7r7l 



2+2^' 



I 2+25 I 2+25' 



L(2e-i)r 



1+^ 



'3+e 



'3+C' 



2^r(^) 



3+i 

1+5 



For any values of K we have 1/2 < < 1. Consequently, the model (|^) develops a massive 
phase for Jm(A) = . The calculation of the bulk free energy leads to the result 



SUSY 
13 



m 



sm 



sin(^)sin(^-; 



f82) 



where we have ^ = K/2, which gives the following expression for the exact vacuum expectation 
value of the perturbing operator 



.2x2,2/ ^ + 2 N ^2+2^ 



TT A 7 ( 



' ^2+2£,' 



2^-1 



X 



2 + 2«'r(i±S)r(|±|)4(i + {)^ 



r2( A)r2(i±S-) sin(^^i^) sinf^ 



•3+£ 



3+5 



'3+5' 



3 + f 



From the exact expression of the vacuum expectation value ( p3D and the definition of the nor- 
malization factor, for more general VEVs of primary fields one obtains: 



-r(^)r(i±|)(i + 



a+ 



^3+5' 



X exp 



T 



-2A 



2 sinh(2(l + ^)t) sinh(2^t) sinh((6 + 2^^) 



^^To describe other perturbations of supersymmetric minimal models, one may have also considered the analytic 
continuation of the model with action (j^). However, for = K/{K + 2) the perturbation corresponds to $1^'^'', 
which is relevant only for K — A. On the other hand, the analytic continuation of action (|^) at the (dual) value 
= {K + 2)/2 with K < A leads to a relevant perturbation, namely ■ For simplicity we will not consider 

this case here. 
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where we introduce the function 



J^{e, i, t) = sinh((l - 9)t){cosh{{e + 3 + 2^^) + cosh((5 + 3)t) - cosh((^ - l)t) 
- cosh((^ - 1 - 2^)t) + (cosh((^ - 3 - 2^)t) - cosh((^ - 7 - 2^^) 
+ cosh((^ + 9 + 2^)t) - cosh((e + 5 + 20t))/2) . (83) 



For K = 2, i.e. ^ = 1 the central charge of the supersymmetric minimal model is csusy = 0. 

(TV. 
13 



Then the perturbing operator reduces to $1^"^^ = I, i.e. the identity operator. In this case, one 



can check that the above relation A-m gives A = — rn^/S which implies fi^^"^ = A, as expected. 

For K = 3, the model with action ( pOD describes an integrable perturbation of the tricritical 
Ising model. On the other hand, the same model can be obtained starting from an integrable 
perturbation of a nonsupersymmetric minimal models, as A^4/5 = Ai^^^. Using notations of 

the previous section, we have the correspondence $12 = $13^"^^- For this case, the VEVs, exact 
relations between the mass scale of the particles and the parameter in the action have been 



proposed in |]rT|. It is straightforward to check that the above results perfectly agree with these 



ones as long as the mass scale m = M of [p 



For K = 6, the central charge of the conformal part becomes c = 5/4, i.e. the resulting 
model Ai^^^ is identified with a critical Zg spin model perturbed by the thermal operator with 
conformal dimension Di = 1/4 denoted e^^-* in In general, Z„ spin models are the natural 
generalizations of the Ising model, when the spin variable takes its values in group Z„ |35|. These 



self-dual models possess critical points which are associated with Zn parafermionic CFTs with 
central charge c = |]35| briefly described in beginning of section 2. Besides the parafermionic 



symmetry, these CFTs possess W(A„_i) symmetry and can be described by the M.n+i{^n-i) 
minimal model at the specific coupling = It has been shown that integrability is preserved 

by adding the perturbing first thermal operator e'-^^ with conformal dimension 2/(n + 2). This 
operator is anti self-dual, i.e. his sign changes under duality transformation. Depending of the 
sign of the parameter characterizing the strength of the perturbation, the perturbed theory will 
be either in ordered or disordered phase. For general values of n, exact relation between this 
parameter and the mass of the particles have been proposed in However, to compare our 
results to the ones associated with the perturbed model Mt^A^), we need to identify the lowest 
kink in both models. Similarly to the case of the coupled models in section 4, for general values 
of ^ here we expect the identification 

m = 2Msin( — (84) 
3 + e 

where M denotes the mass of the lowest kink. Taking ^ = 3 and using the substitution above 



in (^) and (q^); is straightforward to check that previous results coincide perfectly with the 



ones in il3|. 



7 Concluding remarks 

In a first part, we have studied in detail two of the simplest dual representations of integrable 
deformations of affine Toda theories, which corresponds to a critical Ising model coupled with a 
Liouville field theory depending on the coupling b. In the deep UV, the effective central charges 
which characterize the behavior of both models are obtained explicitly. Also, we propose the 
exact relations between the parameters in the Lagrangians (UV data) and the mass scale of 
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the particles (IR data), which are necessary to perform properly the TBA analysis. Various 
vacuum expectation values, namely the one associated with the primary field {e'^'^{x)) and the 
first descendant field (ipipe'^'^lx)) , as well as the bulk free energy are obtained explicitly. All 
previous results are shown to exchange under a weak-strong coupling duality 6 <-> 1/6, confirming 
the duality relation between both Lagrangian representations conjectured in p, ^. 

In the second part, we consider various applications of these results. From the model (H) 
for n = 2, we were able to describe several kind of coupled models. First, we studied in detail 
the case 6^ = —p/2p' or 6^ = —p'/2p from which we obtain the coupled minimal models ( pSf ) or 
(p^, respectively. The main difference between both actions is the presence of an extra term 
in (|53|), associated with the $13 primary operator. It comes from the counterterm in the real 
coupling QFT. In both cases, we propose exact relation between the mass of the lowest breather 
m and the mass of the lowest kink M , as well as the exact relation between the parameters A or 
A and M. For a critical Ising model coupled with unitary minimal models, the QFTs (|53D and 
( ^4|) are found to be massive. Exact VEVs and bulk free energies are obtained. Some special 
cases have been independently checked and shown to agree using different methods^. Among 
these coupled models, we propose an identification between a special case (p = 0) of the SU (8)2- 
HSG model for a finite resonance parameter and a subsector of the Ising- Tricritical Ising models 
coupled through energy-energy. In particular, we find that the lowest breather disappears from 
the spectrum, unstable as expected p3 



Secondly, we propose to study the vicinity of the critical point P in the phase diagram of 
the coupled Ising-XY model depicted in figure 1. Three candidates, integrable perturbations of 
a c = 3/2 CFT, can be considered. The supersymmetric sine-Gordon suggested in or the 
models (|ID or for imaginary coupling and n = 2. Here we focused on the QFT (|76) obtained 
from (0) for n = 2 as the other results can be found in The model is a massive QFT, 



corresponding to a first order phase transition along the line PT in figure 1. Exact VEV of the 
electric/magnetic operator (O±i,o) = (C^cii) is calculated. 

From the model (|I]) for n = 2, we study in detail the NS sector of the QG restriction associ- 
ated with 6^ = —K/K + 2, which leads to the (^f^^^ integrable perturbation of = 1 unitary 
minimal models. We show that results agree with the known ones obtained for the perturbed 



Tricritical Ising model [|TT| for K = 2, and the perturbed parafermionic CFTs |T3[ for K = Q. 
Finally we would like to add a few remarks: 



• Scattering theory of solitons for imaginary coupling. 

For imaginary coupling both models studied here admit a quantum group symmetry based on 
Uq{A^^\{)^2)) and Wg(-B^^^(0, 1)), respectively. Then, it would be interesting to construct the 
corresponding i?-matrices in such a way to obtain the exact S'-matrices for the scattering of soli- 
tons. The identification of the scattering amplitudes for the lowest breathers with the S'-matrix 
for quantum particles in the real coupling case should confirm eq. (^). Also, the restricted R- 
matrix with respect to the different subalgebras should provide the scattering amplitudes in the 
coupled models described here. Understanding the vacuum structure of all the models described 
above would also allow one to fix the form of the prefactor associated with the QG restriction. 



• Homogeneous sine-Gordon 5*^7(3)2 model. 

As we mentioned above, the perturbing field ^per-t — Q{x) with conformal dimension 3/5 can be 
written as ^pf.rt{x) = e{x)^i2{x) + p^i3{x) where p is a c— number. For a finite resonance param- 
eter ||2T|, |2^ the 5'[/(3)2-IISG model contains only two self-conjugate solitons with masses M and 
M'. Then, it would be rather interesting to obtain the exact relation between the ratio M/M' 
and the parameter p. With the A-M relation (M is the lightest mass of the two self-conjugate 
solitons) obtained from eq. (|6D for ^ = 2/3, i.e. |A| = 0.2790872531. ..M^/^ and using (|58|) and 
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it would be straightforward to obtain the expectation values {Oq'1{x)) and {^pert{x)) in the 
HSG-S't/(3)2 model for a finite resonance parameter. 

Furthermore, it is well-known that any correlation function of local fields {O{r)O{0)) in the 
short-distance limit can be reduced down to one-point functions {0'{r)) by successive applica- 



tion of the operator product expansion H9[. Along the line of pO[, using the exact results 



obtained here as well as the three-point functions calculated in [^| and the relation between the 
HSG-S'f/(3)2 and the coupled Ising-tricritical Ising models proposed in section 4, the UV behav- 
ior of the two-point correlation functions can be studied. On the other hand, the long-distance 
behavior (IR) of the two-point correlation functions in the HSG-S'f/(3)2 model has been studied 
in detail in [g^ using the so-called form-factors approach. Then, it would be interesting to com- 
pare both results which should agree in a good approximation in the intermediate region Mr ~ 1. 

• Ramond sector of = 1 superconformal unitary minimal models. 

In this work we mainly focused on the NS sector of both models. However one may be interested 
in the R sector, for instance, the VEVs of the primary fields (cre"'^),-. In this case we have to 
consider the reflection amplitudes in the R sector of the SL field theory instead of the NS one. 
This reflection amplitude has been calculated in with the result 



7CK , 1 + 6^l r(i - (a + Q')b)TC^ - (a + Q')/b) 



W " 2 



Ta + {a + Q')b)Ta + ia + Q')/b) 



where 2Q' = b + 1/6. If we define Gi^\a) = ((xe"^)^ for the VEVs of the R primary fields in 
both models, instead of (^) we have 

Gf\-a) = Rf\-a) Gi''\-2Q' + a) . 

For the QFT (|l|), the minimal solution of this reflection equation together with (|32D can be 
obtained as before. The result can be written in terms of the NS one, i.e. 



(ae'^'^),=i = (a),=iG'fii(a) exp 

with 



dt sinh(a6t) sinh((a6 - 2Qb)t)f{t, b) 
Tsinh(t) sinh(t62) sinh(i72(l + b'^)t)) 



fit, b) = sinh((62 + l)t) - sinh(6^t) - sinh(t) 

and {a)r is the VEVs of the Ising spin field in the QFT The result for the QFT (|^) is 
straightforward. 

Considering specific values of the coupling constant in the QFT (J^) it is then possible to 
obtain the VEVs (o"$rs) in the coupled minimal models studied in sect. 4, which provides the 
leading term of the two-point function between operators which belong to the critical Ising model 
and minimal model, respectively. Also, from the QFT one obtains the primary operators in 
the R sector of the perturbed A^ = 1 superconformal minimal models. 

To conclude, as we mentioned briefly in the Introduction the models studied here belong to 
a more general family of deformed Toda models based on Lie superalgebras. In these cases, 
previous analysis can be performed along the same line. At specific values of the coupling, 
beyond describing different kind of coupled models, these series are identified with integrable 
perturbation of CFTs associated with WB{0,r) algebras with central charge 

CwB(o,r) = (r + 1/2)(1 - 2r(2r - 

pp 
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for p,p' > 2r — 1. We intend to return to these models in the future. 
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Appendix A: Normalization of the VEV of the first descendant field 

To calculate the VEV of the first descendant field {ipip exp(a(y9))T-, we considered H^-^a) defined 
in eq. (|^) as the minimal solution of certain reflection equations. However, we used the exact 
result for the bulk free energy to fix its overall coefficient. It is then rather important to possess 
an independent derivation for this coefficient, based on the "resonance condition" 0. 

For instance, let us consider the analytic continuation h ^ i(3 and /i — > — /i of model (|l]) with 
n = 2. In the free field theory, the composite field ■?/''?/' exp(iQ;y9) is spinless with scale dimension 

D = A + A = l + al (85) 

For generic values of the coupling f3 some divergences arise in the VEVs of the descendant fields 
due to the perturbation in (|l]) with imaginary coupling. They are generally cancelled if we 
add specific counterterms which contain spinless local fields with cutt-off dependent coefficients. 
If the perturbation is relevant, a finite number of lower scale dimension couterterms are then 
sufficient. However, this procedure is regularization scheme dependent, i.e. one can always add 
finite counterterms. For generic values of a this ambiguity in the definition of the renormalized 
expression for the descendant fields can be eliminated by fixing their scale dimensions to be (|85[) . 
Here, this situation arises if two fields, say Oa and Oq/, satisfy the "resonance condition" ||51|| 

D^ = Da' + k{l- (3^) +l{2- /3^) with {kJ)eN (86) 

associated with the ambiguity 

Oa^Oa + t^ii'O^. (87) 

In this specific case one says that the renormalized field Oa has an (fc|/)-th resonance [^] with 
the field Oa'- In particular, the first descendant field ipip e'x.Y){ia(p) has a (1|0) resonance with 
the field exp(i(Q; — (3)(p) at a = 0. If one looks at the short distance behavior in r = |a;| — > of 
the two-point function 

(e^"i'^(x) V?(0)e^°^'^(0)),=i with ai + a2 = a , (88) 

the contributions brought by Hr=i{ia) and {exp{—ij3(p))r=i have the same power behavior in 
r. Furthermore, H.r=i{ia) and the coefficient in front of the second VEV in the short distance 

^^I am very grateful to Al. Zamolodchikov for suggesting this check. 
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expansion of (^) both exhibit a pole at a = 0. By analogy with the method used in [WL 
require that the divergent contributions compensate each other. This leads to the relation 



^ {exp{-i(3ip))r=i . 



we 



^9) 



It is straightforward to check that (0) with r = 1 and b —>■ i(3 indeed satisfies this requirement. 
In particular, for r = 1, 2 this gives a further support to {il)%l)e~^''')r and the exact bulk free 
energy proposed for both models. 



Appendix B: Numerical values for coupled minimal models 

• Two energy-spin coupled Ising models. 

It corresponds to p = 3, = 4 in (|53|) , i.e. ^ = 3/5 in (^). In this case $12 is the spin operator 
of the second model in (|53D with conformal dimension A12 = 1/16 whereas $13 with conformal 
dimension A = 1/2 is the energy operator. We obtain 

($i2(x)) = 1.281110557.. .M^/^; 
{e{x)^i2{x)) = 7.253910604.. .M^/^ 

where the parameter A is related to the mass of the lowest kink by = 0.07660622552... M^/^. 

• Two energy-energy coupled Ising-tricritical Ising models. 

The case p = 4, p' = 5 in (|53D describes a critical Ising model which interacts with a tricritical 
Ising model through their leading energy density operators. For the second model $12 has 
conformal dimension A12 = 1/10. It corresponds to (3'^ = 2/5 i.e. ^ = 2/3 in (0). This model 
also contains the sub-leading energy density operator $13 with A13 = 3/5 ("vacancy operator"), 
two magnetic operators $22 with A22 = 3/80, $21 with A21 = 7/16 and $14. We have for instance 

1.495279412.. .M^/^; 
1. 133076821. ..M^/^°; 
4.478886063.. .M^/^ 

where the parameter A is related to the mass of the lowest kink by A^ = 0.07788969384.. .M^/^ 



• Two Ising-As RSOS coupled models. 

It corresponds to the choice p = 5, p' = 6 in (0). The RSOS model possesses a primary 
operator $12 with conformal dimension A12 = 1/8 and $22 with A22 = 1/40. Taking = 5/7 in 
(p8D, we obtain 

($22(0;)) = 1. 081616064.. 
($i2(x)) = 1. 673156742.. .M^/^; 
(e(x)<l>i2(x)) = 3.478873566.. .M^/^ 

where the parameter A is related to the mass of the lowest kink by A^ = 0.07848322325.. .M3/2. 
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• Two energy-spin coupled Ising-tricritical Ising models. 

For p = 4, p' = 5 in (|5^) , the Ising energy operator is coupled to the subleading spin operator of 
the tricritical Ising model. For ^ = 2/3 in (|67D we have for instance ($22(3^)) = 1.2148...M^/^° 
where the parameter A is related to the mass of the lowest kink by = 0.0319842.. .M^/^ 



• Two energy-energy coupled Ising-3-state Potts models. 

The case p = 5, p' = 6 in (^) describes a critical Ising model coupled to a critical 3-state Potts 
models through their energy density operator e and $21 with conformal dimension A21 = 2/5, 
respectively. It corresponds to ^ = 5/7 in (|67D . The 3-state Potts model also contains 
the primary operator $23 - the spin operator - with A23 = 1/15. We obtain for instance 
($23(2;)) = 1.3378...M^/^^ where the parameter A is related to the mass of the lowest kink 
by A^ = 0.0420507.. 
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